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Formulas of Drought Indices 


1- Meteorological drought 
A: Rain Based-drought indices (Salehnia et al., 2017): 


e SPI (Standardized Precipitation Index, McKee et al. 1993, 1995 
e DI (Deciles Index), Gibbs and Maher, 1967 

e PN (Percent of Normal Index), Willeke et al. (1994) 

e CZI (China-Z Index), Wu et al. (2001) 

e MCZI (Modified CZI), Wu et al. (2001) 

e EDI (Effective drought Index), Byun and Wilhite (1999) 

e RAI (Rainfall Anomaly Index), van Rooy (1965) 

e ZSI (Z-score Index), Palmer (1965) 


B: Other meteorological drought indices: 


e PDSI (Palmer Drought Severity Index), Palmer (1965) and Dehghan et al., 
2020 

e PHDI (Palmer Hydrological Drought Index), Palmer (1965) 

e SPEI (Standardized Precipitation Evapotranspiration Index), Vicente-Ser- 


rano et al., 2010 and Salehnia et al., 2020 


e RDI (Reconnaissance Drought Index), Tsakiris and Vangelis, 2005. 


2- Agricultural drought indices 


e ARI (Agricultural Rainfall Index), Nieuwolt, 1981 
e SMDI (Soil Moisture Deficit Index), Narasimhan and Srinivasan, 2005 
e ETDI (Evapotranspiration Deficit Index), Narasimhan and Srinivasan, 


2005 
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3- Hydrological drought indices 


1.1 


e SWSI (Surface Water Supply Index), Garen, 1993 
e SDI (Streamflow Drought Index), Nalbantis and Tsakiris, 2009 


SPI (Standardized Precipitation Index) 


The SPI stands as the predominant drought index on a global scale, 
utilized for monitoring and comprehensive analysis, as highlighted 
by Karabulut (2015). Widely acknowledged, the SPI serves as a cru- 
cial instrument for delineating meteorological drought patterns, as 
evidenced by the works of Hayes et al. (1999) and Deo (2011). The 
pioneering work of McKee et al. (1993, 1995) encompassed the es- 
tablishment of SPI across multiple temporal scales (1, 3, 6, 12, 24, 
and 48 months). The resultant SPI values span from +2.0 to -2.0. 
Given the possibility of precipitation data conforming to a gamma 
distribution, the calculation of SPI entails the utilization of the prob- 


ability density function inherent to the gamma distribution: 
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1.2 


g(x) =~ x*1eB for x >0 


In this context, the gamma function is denoted by T, where x repre- 
sents the quantity of precipitation (x >0), a (> 0) stands as the shape 
parameter, and B (> 0) represents the scale parameter. For further 
insights, refer to the comprehensive studies conducted by Edwards 


and McKee (1997) as well as Dogan et al. (2012). 


PN (Percent of Normal Precipitation Index) 


The PN index, delineated by Willeke et al. (1994), represents the 
percentage of normal precipitation attributed to a particular loca- 
tion. This metric, PN, holds the versatility to be computed over var- 
ious temporal scales such as monthly, seasonal, and annual periods. 
Notably, when applied to a singular geographic region or specific 
season, PN has demonstrated its efficacy in elucidating drought 
conditions, as demonstrated by the research of Hayes (2003). The 


calculation of PN is as follows: 
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PN ="*/, x 100 
Here, Pi denotes the precipitation during the time increment la- 
beled as "i," while P represents the mean normal precipitation for 


the duration of the study period. 


DI (Deciles Index) 
The deciles index (DI) was defined by Gibbs and Maher (1967). 


Monthly historical precipitation data are sorted from lowest to 
highest and divided into ten equal categories or deciles, so a given 
month’s precipitation is placed into historical context by decile. 

To compute the deciles index in DMAP, we first calculate the deciles 
for each month and subsequently identify any specific month-year 


in which the deciles meet certain criteria. 


DEC, = Per(Pn,N) 
Dli; =n if DECr-1 < Pij Ss DEC, 


The term DEC, represents the deciles associated with the value of n, 
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where n takes on values such as 10, 20, ..., 90. Here, Per signifies the 
percentile function. The Pm corresponds to the rainfall observed in 


month m. 


1.4 EDI (Effective Drought Index) 
The Effective Drought Index (EDI) is computed at a daily time step 


and follows a standardization process akin to SPI. Byun and Wilhite 
(1999) pioneered the development of EDI, aiming to address certain 
shortcomings observed in other indices like SPI. The EDI's numerical 
range spans from -2.5 to 2.5. When the EDI falls within the range of 
-1 to 1, it signifies near-normal conditions, whereas values below - 
2 indicate severe drought. The calculation of EDI is accomplished 


through the application of the subsequent equation: 


—————— 
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ie 5 [or Fem) 


DS 


DEP, = EP; — MEP 


_ DEP, 
EDI; =" '/sp (DEP) 


Where DS represents the count of days used for summing precipi- 
tation in the calculation of drought severity (usually 365 or 15), SD 
stands for standard deviation, and MEP stands for the mean of EP. 

The DMAP calculates EDI on daily scale then you can use Monthly, 


Yearly, or seasonally to convert daily EDI to this time scales. 


1.5 CZI (China-Z Index) 
In 1995, the National Climate Center of China introduced the CZI as 


a substitute for the SPI (Ju et al., 1997) specifically designed for sit- 
uations where mean precipitation adheres to the Pearson Type III 


distribution. The CZI is computed as follows: 


6 Csi 6 C3 
Sl Sl 
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aay — X1)° 


n xop 


Csi = 


where i is the time scale of interest, j is the current month, Csi is the 
coefficient of skewness and qj; is the standardized variation. Fur- 


ther details can be found in Wu et al., 2001. 


1.6 MCZI (Modified CZI) 
The adjusted MCZI employs the aforementioned formula but re- 


places the mean precipitation with the median precipitation. 


1.7 ZSI (Z-Score Index) 
The Z-Score is sometimes mistaken for SPI. However, it bears a 


closer resemblance to CZI, although it doesn't necessitate fitting 
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precipitation data into either Gamma or Pearson Type III distribu- 
tions: 


P; — P 
SD 


ZSI = 


Here, P represents the average monthly precipitation (mm), Pi 
stands for the precipitation in a particular month, and SD denotes 


the standard deviation across any time scale. 


RAI (Rainfall Anomaly Index) 

The RAI involves two phases: positive anomalies and negative 
anomalies. Initially, precipitation data are sorted in descending or- 
der. The top ten and bottom ten data points are then averaged 
separately to create these phases. The positive and negative 


anomalies are computed using equations 6 and 7, respectively: 


B (p — p) 
RAI =3 x Gop) 
_—_ a eae) 
RAI = —3 x Gs 


In this context, RAI stands for Rainfall Anomaly Index. Here, p 
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represents the actual precipitation (mm) for each year, p is the 
long-term average precipitation, and M is the mean of the top ten 
values of p for the positive anomaly, as well as the mean of the 
bottom ten values of p for the negative anomaly. 

Reference for 1.1 up to 1.8: Estimation of meteorological drought indices 


based on AgMERRA precipitation data and station-observed precipitation 


data 


1.9 KBDI(Keetch-Byram Drought Index) 


Calculating the index requires data for daily maximum tempera- 
ture and precipitation. Following that, the KBDI was computed us- 


ing the subsequent methodology: 


[800 — KBDI,_,] x [0.968 x Exp(0.875T max + 1.5552) — 8.30] 
= EET, x 
1+ 10.88 x Exp(—0.001736R) 


DF 1073 


KBDI,_4 if P=0 and Tmax < 6.78C 
KBDI,_, + DF if P=0 and Tmax > 6.78 C 


KBDI,= 4 KBDI,., +DF if P>0O and ` P < 5.1mm 


KBDI; + DF if P>0 and X P> SA mm 
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KBDI, = KBDI,_, — 39.37 X X} p 


The drought factor (DF) on a given day in the metric system, as pre- 
sented by Janis et al. in 2002, is calculated based on several factors. 
Specifically, it depends on the daily maximum temperature (Tmax:) 
in degrees Celsius, the average annual rainfall (R) in centimeters for 
the region, the Keetch-Byram drought index for the previous day 
(KBDhk-1), and the daily precipitation (P+) in millimeters (Janis et al. 
2002). The calculation involves intricate interplays between these 
variables. Notably, the initialization of the Keetch-Byram drought 
index (KBDI) is of paramount importance. Traditionally, this initiali- 
zation occurs after a series of consecutive rainy days when soil sat- 
uration is attained. However, it's important to acknowledge that the 
field capacity of arable land, which signifies its ability to absorb wa- 
ter, can differ significantly due to variations in soil composition and 


changes in vegetation. 
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Referenc: Predictive value of Keetch-Byram Drought Index for cereal yields 


in a semi-arid environment 


1.10 PDSI 


1.11 PHDI 


Calculating these two indices involves a substantial and intricate 
process. It's advised to consult the original paper (Meteorological 


drought Palmer 1965 Washington DC ) for detailed information. 


Additionally, our calculation of the indices is directly based on this 


paper: A tool for calculating the Palmer drought indices 


1.12 SPEI 


In the computation of SPEI, we employed both the classical approx- 
imation by Abramowitz and Stegun (1965) and the method pro- 
posed by Vicente-Serrano et al. (2009). The SPEI, accumulated over 
various time scales, is founded upon a monthly climatic water bal- 


ance derived from the difference between precipitation and 
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Evapotranspiration (ETO). This balance is then refined through the 
application of a three-parameter log-logistic distribution. The accu- 
mulation of deficit or surplus in the climatic water balance (D) at 
different time scales is ascertained by evaluating the dissimilarity 
between daily precipitation (P) and ETO for a given day "i": 
Dj = P; - ETo; 


The probability density function of a three-parameter log-logistic 


distributed variable is expressed as below: 


poy EED xfs EY] 


Here, a, B, and y represent the scale, shape, and origin parameters, 


respectively, for D values within the range (y > D < °°). The parame- 
ters for the Pearson Ill distribution can be acquired using the fol- 


lowing method: 


2W1 — Wo 


p= 


6W, — Wo — 6w 
1 0 2 


_ (Wo — 2w,)B 
pee 
ra + /pra z 1/9) 
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p= m E)E) 


Here, [(B) denotes the gamma function of B. The probability distri- 


bution function of the series "D" is expressed as follows: 


In the final stage, by utilizing the value of F(x), the SPEI can be de- 
termined as the standardized value of F(x). The calculation of the 
SPEI equation is as follows: 


Co + CGW + CW? 
1+d,W+d,W2+d,W? 


SPEI = W 
Here, for cases where P < 0.5, W is defined as the square root of -2 
times the natural logarithm of P, where P represents the probabil- 
ity of surpassing a specific value of "D." The constants involved are 


as follows: CO = 2.515517, C1 = 0.802853, C2 = 0.010328, d1 = 


1.432788, d2 = 0.189269, and d3 = 0.001308. 


Reference: Rainfed wheat (Triticum aestivum L.) yield prediction 
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using economical, meteorological, and drought indicators through 


pooled panel data and statistical downscaling 


RDI 
RDI (Reconnaissance Drought Index), developed by Tsakiris and 
Vangelis in 2005, compares two cumulative measures: precipita- 
tion and potential evapotranspiration. For a given month "m", this 
calculation for a yearly index is represented by the following equa- 
tion: 
D 
Dm = nA 
i=Jan i 


The Reconnaissance Drought Index (RDI) is characterized by two 
variations: the normalized RDI and the standardized RDI. The nor- 
malized RDI serves as the default option within DMAP. 

The Normalized Reconnaissance Drought Index (RDI) is computed 


using the following equation: 


Dm 
RDI =—-1 
m 
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The D is average of D for all years. The Standardized Reconnais- 
sance Drought Index (RDI) is computed using the following equa- 


tion: 


m 


RDI m = 


The D is average of D for all years and ø is standard deviation of 
D for all years. 

Reference: Regional Drought Assessment Based on the Reconnais- 
sance Drought Index (RDI 

SMDI 

For the computation of SMDI, it is necessary to employ the long- 
term median, maximum, and minimum monthly (or seasonal or 
yearly) soil water values. The following formulas are utilized for 


this purpose: 


SWij -MSW 100 if SwW.. < MSW 
MSW; — mins W, ia aes 
SD, ; = 
SW; ; — MSW, 


—~ xX 100 if SW,;; > MSW; 
maxSW; — MSW; i J 


SD; 
SMDI; = 0.5 x SMDIj- += 
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The SMDI values will range from -4 to 4 for comparison with PDSI 
values. 

Here, MWS; represents the long-term median of water stress for 
month j, maxMWS; stands for the long-term maximum water 
stress of month j, minWS; signifies the long-term minimum water 
stress of month j, and WS corresponds to the monthly water 
stress. The subscripts i and j are employed for years and months 
respectively. The water stress anomaly varies from -100 to +100 
for each month, denoting exceedingly dry to exceptionally wet 


conditions. 


Reference: Development and evaluation of Soil Moisture Deficit 
Index (SMDI) and Evapotranspiration Deficit Index (ETDI) for agri- 


cultural drought monitoring 


ARI 
ARI serves as a dependable indicator of monthly water equilib- 
rium, specifically designed for assessing agricultural droughts. The 


calculation of ARI follows this methodology: 
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ARI = 2 
~ PET 


Here, P represents the monthly precipitation (mm/month), and 
ETO signifies the reference evapotranspiration (mm/month). The 
presence of drought within a specific month is denoted by ARI val- 
ues of 40. A range of ARI values between 40 and 200 indicates fa- 
vorable conditions for vegetation growth and agricultural produc- 
tivity. An ARI exceeding 200 for a month signifies a period of ele- 
vated moisture, indicating wet conditions. 

Reference: Understanding Dry and Wet Conditions in the Vietnam- 


ese Mekong Delta Using Multiple Drought Indices: A Case Study in 


Ca Mau Province 

ETDI 

The ETDI is computed similarly to the SMDI, yet it relies on the de- 
viation of water stress from its long-term mean. Here, the monthly 
water stress is established through the interplay of potential and 


actual evapotranspiration. 


_ PET — AET 


WS 
PET 
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Wij 00 if WS,, < MWS 
MWS, — minWS; AE 
WSA; j = 
WS; ; — MWS; 


I X100 if WS, > MWS; 
maxWS; — MWS; WS j 


WSA; 


ETDI; = 0.5 x ETDIj_, + 0. 
The ETDI values will range from -4 to 4 for comparison with PDSI 


values. 


Reference: Development and evaluation of Soil Moisture Deficit 
Index (SMDI) and Evapotranspiration Deficit Index (ETDI) for agri- 


cultural drought monitoring 


1.17 SWSI 
A prevalent hydrological drought index consists of a series of prac- 
tical individual measures, known as the Surface Water Drought Se- 
verity Index (SWSI). The formula used to compute this index, as 


per Garen's investigations, can be described as follows: 


100 x P;;— 50 
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P,j = GammaCDF (StreamF low; j) 
The GammaCDF represents the cumulative distribution function of 
the gamma distribution fitted onto month j across all years. In this 
context, StreamFlow;; denotes the streamflow during month j of 
the year i. 
Reference: Comparing surface water supply index and streamflow 
drought index for hydrological drought analysis in Ethiopia 
SDI 
The Streamflow Drought Index can be computed using the follow- 
ing formula when applied to monthly data. However, if you pro- 
vide daily data, the tool will automatically aggregate it into 


monthly values before performing the calculation. 
3K 
Vik = X SF; i = 1,2,..,endYear j=1,2,..,12 k 
j=1 


= 1,2,3,4 


Vik — Ve 


SDIix = Eo 


The premise is that a time series of monthly streamflow volumes 


is available, denoted as SFij. Here, i indicates the hydrological year, 
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and j represents the month within that hydrological year (where j 
= 1 corresponds to October and j = 12 corresponds to September). 
In this context, V; and SDx stand for the mean and standard devi- 
ation of cumulative streamflow volumes observed during the ref- 
erence period k, established over a substantial time span. Within 
the framework of DMAP, the resulting output will encompass 
these four categories: Oct-Dec, Oct-Mar, Oct-Jun, and Oct-Sep 


(Yearly). 


Reference: Assessment of Hydrological Drought Revisited 


e 
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